Best 

Available 

Copy 


U.D.C.  621.396.969.571.11  :  521.5  :  629. 19. 096  :  518.5 


ROYAL  AIRCRAFT  ESTABLISHMENT 
Technical  Report  67199 
JUly  1967 

ON  THE  USE  OF  RANCHI  AND  RANGE  RATE  MEASUREMENTS 
IN  THE  NEW  R.A.E.  ORBIT  DETERMINATION  PROGRAM 

by 

R.  H.  Gooding 


SUMMARY 

The  new  orbit  determination  program  PROP  can  use  range  and  range  rate 
observations  as  well  as  observations  of  direction.  The  orbit  of  Ogo  2 
(1965-81A)  has  been  determined  at  epoch  1965  November  10.0,  using  Minitrack 
observations  and  Goddard  Range  and  Range  Rate  observations.  The  results 
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show  reasonable  agreement  between  observations  from  the  two  systems  and 
suggest  that  random  errors  in  the  range  and  range  rate  observations  are  of 
the  order  claimed. 

This  paper  was  presented  at  the  8th  COSPAR  International  Space 
Science  Symposium,  London,  July  1967. 
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1  INTRODUCTION 

Past  determinations  of  earth-eatellite  orbits  at  the  Royal  Aircraft 
Establishment  have  been  restricted  to  the  use  of  directional  observations 
only,  i.e.  to  observations  of  right  ascension  and  declination,  azimuth  and 
elevation,  or  a  pair  of  direction  cosines .  During  the  last  two  years,  however 
a  new  computer  program  has  been  developed  which  permits  the  use  of  other  types 
of  observation,  including  range,  range  rate  and  the  rates  of  change  of 
directional  data. 

The  first  application  of  the  new  program  has  been  to  the  orbit  of 
Ogo  2  (1965-6IA)  over  a  period  of  four  days  in  November  19^5.  This  satellite 
was  tracked  by  two  RASA  systems  of  stations  -  the  Minitrack  network  and  the 
Goddard  Range  and  Range  Rate  System.  Thus  an  opportunity  was  available  to 
test  the  ability  of  the  new  program  to  deal  with  observations  of  p  and  p 
(range  and  range  rate)  and  to  see  how  the  accuracy  of  orbital  parameters 
determined  from  Goddard  p  and  p  observations  compared  with  their 
accuracy  when  determined  from  Mini  track  observations.  There  was  also  the 
opportunity  to  form  some  assessment  of  the  accuracy  of  the  p  and  p 
observations  themselves. 

2  THE  GODDARD  RANGE  ARP  RANGE  RATE  SYSTEM 

The  GRARR  System  was  developed  to  supplement  the  Mini  track  network, 
which  consists  of  about  a  dozen  interferometer  stations.  Minitrack  observa¬ 
tions  are  very  suitable  for  close-earth  satellites,  but  for  satellites  at 
greater  distances  they  suffer  from  a  weakness  which  is  inherent  in  all 
directional  observations:  viz.  that  a  fixed  error  in  angle  generates  an 
error  cone  and  hence  is  equivalent  to  a  positional  error  which  is  proportional 
to  range. 

The  GRARR  System  consists  of  stations  at  Carnarvon  (Australia), 

Rosman  (North  Carolina),  Tananarive  (Malagasy),  Fairbanks  (Alaska)  and 
Santiago  (Chile).  The  first  three  were  constructed  first,  and  data  from  the 
other  two  have  not  been  received.  Each  station  operates  at  either  VHF  or 
S-Band  and  is  capable  of  measuring  range  and/or  range  rate.  (Stations  can 
also  measure  direction,  but  errors  may  exceed  0.1*,  so  that  this  is  not 
normally  very  useful.)  The  Ogo  2  observations  received  from  NASA  used  the 
S-Band  mode  of  operation,  which  is  the  more  accurate  and  is  likely  to  be  the 
more  extensively  u?>ed.  The  VHF  antenna  is  required  even  in  the  S-Band  mode, 
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however,  since  It  Is  responsible  for  Initial  acquisition,  by  locking  onto 
the  136  MHz  Mini  track  beacon  In  the  satellite.  The  S-Band  antenna  is 
initially  slaved  to  the  VHF  antenna  but,  after  acquisition  by  the  S-Band 
receiver,  it  adopts  an  autotrack  mode  of  operation. 

Range  is  measured  by  use  of  a  ground-transmitted  100  kHz  side  tone  fre¬ 
quency  and  a  series  of  lower  frequency  tones  to  resolve  ambiguities.  Range 
rate  is  measured  by  counting  cycles  of  the  two-way  doppler  shift. 

The  performance  objectives  for  the  GRARR  System  specify*  S-Band 
accuracies  (rms  tolerances)  of  15  m  in  p  and  0.1  b/b  in 

3  000  2  OBSERVATIONS 

Minitrack  observations  on  punched  cards  have  been  received  (at  R.A.E.) 
for  the  period  November  9-11  1965,  and  ORARR  observations,  also  on  punched 
cards,  for  the  period  November  7-12  1965.  It  was  decided  to  work  with  about 
60  of  each  type  of  observation,  covering  the  period  November  8-11,  and  to 
refer  orbital  parameters  to  a  single  epoch,  namely  19&5  November  10.0. 

The  observations  were  weighted  in  the  orbit  determination  program, 

using  the  following  a  priori  estimates  of  error  (s.d.)t  for  the  Minitrack 

observations,  0.00029  in  the  direction  cosines,  equivalent  to  about  1 1  in 
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angular  measure,  as  suggested  by  previous  R.A.E.  experience  ;  for  the  GRARR 
observations,  150  m  in  p  and  ^  m/s  in  p,  assessing  these  much  more 
conservatively  than  the  performance  specification. 

4  ORBITAL  MODEL 

The  model  for  the  new  R.A.E.  orbit  determination  program  (known  as  PROP) 
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has  been  described  by  Merscn  .  It  takes  account  of  short-periodic,  secular  and 
long-periodic  perturbations  due  to  the  zonal  harmonics  of  the  earth*s  gravita¬ 
tional  field  and  to  atmospheric  drag. 

There  is  sane  flexibility  in  the  choice  of  the  particular  set  of 
orbital  parameters  to  be  fitted  by  PROP.  The  following  set  was  chosen: 
e  (eccentricity),  i  (inclination),  fi  (right  ascension  of  the  node), 

«*>  (argument  of  perigee),  Mq  (mean  anomaly),  M1  and  these  are  as  sinned  to 
be  epoch  values  in  each  case.  The  last  two  parameters  are  such  that  the  mean 
anomaly  at  time  t  (relative  to  epoch)  is  given  -  apart  from  short-periodic 
and  long-periodic  perturbations  -  by 

M=MQ  +  M1t4-M2t2 
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Thus  Hj  is  the  neac  motion  at  epoch;  it  is  related  to  a,  the  aemi-maj or 
axis,  by 

a  -  (n/M^)1/5  -  i  J2  R2  Gi/M^)”1/5  (1  -  l£  Bin2  DO  -  e2)'1* 

Here  p,  J2  and  R  are  familiar  earth  constants,  for  which  the  following 
values  have  been  taken: 

p  -  398602  krn^/s2,  J2  -  TO"6  X  1082.68,  R  -  6378.163  km  . 
5  RESULTS 

The  determination  of  orbital  parameters  at  the  selected  epoch  was 
carried  out  three  times:  first,  using  39  Minitrack  observations;  second, 
using  27  range  and  28  range  rate  observations;  finally,  using  all  11 4 
observations  together.  The  results  are  summarized  in  Table  1.  It  is  to  be 
remembered  that  there  are  only  7  independent  parameters;  the  semi -major 
axis  a  is  a  derived  parameter. 

Table  1 


Case  1 

Case  2 

Case  3 

a  (km) 

7344.8382  ±0.0008 

7344.8427  ±0.0002 

7344.8414  ±0.0002 

e 

0.075948  ±4  x  10"° 

0.075950  i£  X  10-t> 

O.075947  ±2  x  10_t> 

i  (deg) 

87.3574  ±0.0007 

87.3591  ±0.0005 

87.3571  ±0.0005 

0  (deg) 

-84.5244  ±0.0006 

-84.5241  ±0.0002 

-84.5245  ±0.0003 

“  (deg) 

92.5875  ±0.0052 

92.6195  ±0.0018 

92.6105  ±0.0024 

MQ  (deg) 

-73.6895  ±0.0048 

-73.7187  ±0.0018 

-73.7110  ±0.0022 

M.j  (deg/d) 

4966.6788  ±0.0008 

4966.6742  ±0.0002 

4966.67.56  ±0.0002 

M_  (deg/d2) 

0.0022  ±0.0006 

0.0024  ±0.0002 

0.0041  ±0.0002 

6  out  of  59 

4  out  of  55 

11  out  of  114 

1.28 

O.56 

1 .06 

Some  explanation  is  necessary  of  the  last  two  rows  of  Table  1 . 

The  quantity  e,  known  as  the  standard  deviation  of  an  observation  of 
unit  weight,  is  a  measure  of  the  goodness  of  fit.  It  is  given  by 
c  a  Vfc/D^,  where  1  is  the  sum  of  the  squares  of  the  final  weighted 
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residuals  and  Df  is  the  number  of  degrees  of  freedom.  If  the  accuracies  of 
the  observations  have  been  correctly  assessed,  and  if  the  errors  in  the 
observations  are  random  and  normally  distributed,  then  the  expected  value 
of  e  is  1 .0.  The  standard  deviations  given  for  the  orbital  parameters  all 
include  the  appropriate  e  as  a  factor;  this  is  why  sons  of  the  standard 
deviations  for  the  third  orbital  determination,  for  which  all  the  observa¬ 
tions  were  included,  are  greater  than  the  corresponding  values  for  the  second 
determination,  for  which  only  about  half  the  observations  were  used  but  for 
which  6  was  only  0.5 6. 

The  computer  program  has  an  automatic  facility  for  rejecting  any  observa¬ 
tion  with  a  weighted  residual  which  exceeds  a  certain  quantity;  in  the  three 
orbit  determinations  here  the  quantity  concerned  is  3®  •  In  the  second 
determination  the  four  observations  rejected  comprised  two  of  range  and  two 
of  range  rate.  Since  in  both  the  first  two  determinations  there  were  several 
observations  with  weighted  residuals  in  the  vicinity  of  $5e  it  is  not 
surprising  that  the  observations  rejected  in  the  third  determination  were  not 
quite  the  sane  as  those  rejected  in  the  earlier  determinations  -  they  comprised 
nine  Minitrack,  two  range  rate  and  no  range  observations.  It  was  because  the 
accepted  observations  were  not  Just  the  sum  of  the  accepted  observations  from 
the  first  two  determinations  that  the  orbital  parameters  do  not  always  lie 
between  the  values  given  by  the  first  two  determinations . 

6  DISCUSS IOH 

An  obvious  conclusion  from  the  values  of  e  in  Table  1  is  that  the 
observations  were  wrongly  weighted.  The  Minitrack  observations  appear  some¬ 
what  less  accurate  than  was  assumed,  while  the  GRARR  observations  are  more 
accurate.  Inspection  of  the  residuals  from  the  second  orbit  determination 
shows  two  things:  first,  that  the  average  residual  in  p  was  only  about 
30  m,  while  that  in  p  was  about  0.3  a^s;  second,  that  the  random  components 
were  actually  much  less  than  this,  a  strong  bias  being  present. 

Biases  of  this  order  cannot  be  regarded  as  very  serious.  There  are 
several  possible  explanations  of  them.  The  accuracy  of  the  station 
coordinates  is  not  known  but  it  is  unlikely  to  be  better  than  30  m;  the 
position  given  for  Rosman  in  fact  differs  by  about  £  km  from  that  for  an 
earlier  survey.  Errors  in  the  orbital  model,  due  in  particular  to  the  neglect 
of  all  but  the  dominant  tea  serai  harmonic  perturbation,  will  certainly  exceed 
30  m;  over  one  pass  of  the  satellite  from  any  station  these  errors  will 
appear  as  bias  in  the  observations. 
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It  is  useful  to  consider  the  differences  in  the  orbital  parameters 
given  by  the  two  Independent  determinations.  These  differences  may  be 
eaqpressed  in  terms  of  the  standard  deviations  of  the  parameters,  where  \m 
choose  the  first  set  of  standard  deviations  since  these  are  larger.  We 
have 


6a  ■ 

0.0045  km  - 

6e  ■ 

2  X  10"6  a 

6i  - 

0*.0017 

2 1®, 

6o  . 

0°.0003  a 

\  0 

6U  a 

0\0520 

6  c  , 

6m  - 

0 

0*,0292  - 

6  a 

6r,  ■ 

0.0046#/d  = 

5|cr 

and  6*t> 

-  0.0002  7 d2 

- 

The  existence  of  and 6  a  differences  again  indicates  bias.  For  semi¬ 

major  axis  the  bias  may  be  due  to  error  in  the  adopted  value  of  |i.  In 

the  case  of  «  and  M  it  must  be  remarked  that  these  parameters  are  very 

o 

highly  correlated,  since  the  position  of  perigee  is  not  well-determined. 

7  CONCLUSIONS 

The  new  R.A.E.  computer  program  can  deal  satisfactorily  with  observa¬ 
tions  of  range  and  range  rate. 

NASA  observations  from  the  Goddard  Range  and  Range  Rate  System  appear 
to  have  accuracies  consistent  with  the  performance  objectives 
(15  m  and  0.1  m/s),  as  far  as  random  errors  are  concerned.  For  observations 
of  Ogo  2  analysed  by  the  R.A.E.  program,  however,  the  residuals  are  of  the 
order  30  m  and  0.3  m/s.  The  effective  bias  in  these  residuals  can  be 
attributed  to  limitations  in  the  orbital  model  of  the  program;  whether  there 
are  also  actual  biases  in  the  observations  cannot  be  decided. 

For  an  orbit  at  an  average  height  of  about  1000  km  orbital  parameters 
4  obtained  from  the  Minitrack  network  and  parameters  obtained  from  the  GRARR 

System  are  of  comparable  accuracy.  In  a  single  conparison  for  the  satellite 
Ogo  2  standard  deviations  for  orbital  parameters  determined  from  Mlnitrack 
observations  were  only  about  half  those  for  parameters  determined  from 
GRARR  observations.  For  orbits  at  greater  heights  it  is  to  be  expected  that 
the  advantage  of  the  GRARR  System  would  become  quite  marked. 
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